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Theory of adiabatic Hexaamminecobalt-Self-Exchange 

R. G. Endres, M. X. LaBute, D. L. Cox 
Department of Physics, University of California, Davis, CA 95616 
(Dated: February 2, 2008) 

We have reexamined the thermally induced Co{N H-ifg^^^^ [Co(II/III)] redox reaction using the 
first principles density- functional-theory method, semiclassical Marcus theory, and known charge 
transfer parameters. We confirm a previously suggested mechanism involving excited state i^Eg) of 
Co(II) which becomes lower than the ground state {*'Tig) in the transition state region. This lowers 
the transition state barrier considerably by ~ 6.9kcxil/mol and leads to a spin-allowed and adia- 
batic electron exchange process. Our calculations are consistent with previous experimental results 
regarding the spin-excitation energy {^Tig) of Co(III), and the fact that an optical absorption peak 
(^Eg) of the Co(II) species could not be found experimentally. Our rate is of order 6 • 10~^(Ms)~^ 
and hence 2 orders of magnitude faster than determined previously by experiments. 



I. INTRODUCTION 

The experimental rate determination and theoretical 
understanding of the Co(A^iJ3)6^+/'^+ redox reaction in 
aqueous solution 



Co{NH3)l+ + Co{NH3)l+ 



Co{NH; 



3^6 



Co{NH^)l+ 



(1) 



has been a great intellectual challenge for several decades. 
Despite substantial effort, the mechanism of the rate is 
still an unsolved problem. Is the reaction spin-forbidden 
and diabatic or are spin-excited states thermally accessi- 
ble which would possibly lead to a spin-allowed, adiabatic 
reaction? 

Experimental studies in the early 1960s suggested an 
extremely slow rate of fc = 1.6- 10-i°(Ms)-i y at 64.5 C 
and 1 M ionic strength using radiocobalt as a tracer. 
Unfortunately, side reactions involving hydrolysis of the 
complexes, which contribute to the rate, have not cor- 
rectly been taken into account. In the 1980s a detailed 
analysis of previous data using the Marcus correlation [p| 
led to a much higher estimate of 1 • 10~^(Ms)~^ [pT, 
and subsequent experiments labeling the ammonia lig- 
ands with ^^N (which can be assayed by NMR of the 
ammine protons and make it possible to trace side reac- 
tions) gave about 6- 10~^(Ms)~i at 40 C and 2.5 M ionic 
strength Q . An estimate of the rate at normal conditions 
25 C and IM ionic strength) gives < 5 • 10""^ (Ms) 
" I, 0, §. These results can be expected to be rather 
trustworthy, because they are, as expected, similar to 
the rate of a closely related and well understood system. 



Coen 



3+/2+ 



which has a rate of 5.2- 7.7- 10-5(Ms)-i §. 
Furthermore, an early but less influential measurement of 
Co{NH3)g^^^^ in liquid NH^ (and hence not suffering 
from hydrolysic side reactions) indicated a similar rate of 
order 10-^{Ms)-^^. 

Briefly, previous results of theory mainly by Buhks 
et al (1978) Q and Newton (1986,1991) |, || indi- 
cate that the reaction involves the ground state species 
low-spin Co{NH^)l+, ^Aig [Co(III,S=0)], and high- 



states, ^Tig [Co(III,S=l] and "^Eg [Co(II,S=l/2], were to 
high in energy to be thermally populated. This led to 
a spin-forbidden (only possible by weak spin-orbit cou- 
pling), diabatic reaction with a rate constant of about 4 
orders of magnitude too small compared to experiment. 
However, work by Larsson et al (1985) stressed that the 
ground state-excited state energy separations at the tran- 
sition state have to be considered, not at the equilib- 
rium geometries. Including Jahn- Teller (JT) stabiliza- 
tion energies for the excited states, this led in particular 
to a substantial energetic lowering of Co(II,S=l/2) rela- 
tive to Co(II,S=3/2). Although this made the reaction 
spin-allowed, their work suffered from inconsistent data 
from different electronic structure codes. Furthermore, 
the rate constant was not determined. 

In this paper we show that the thermal activation bar- 
rier for the spin-allowed reaction between Co(III,S=0) 
and Co(II,S=l/2) [total spin S in units of ?i] is drastically 
lowered (~ 6.9kcal/mol) compared to the spin-forbidden 
reaction between Co(III,S=0) and Co(II,S=3/2). Un- 
like previously believed, the JT-distorted spin-excited 
Co(II,S=l/2) and the ground state Co(II,S=3/2) are 
near degenerate. Hence the reaction is spin-allowed and 
adiabatic. Our results are based on the density functional 
theory (DFT)Q code SIESTA|14) which we use to cal- 
culate spin-excited states of isocharge molecular species 
and potential energy surfaces (PESs). Since DFT is a 
ground state theory, excitation energies obtained from 
ground state energy differences of species with different 
total spin are expected to be rather good. The DFT 
method has the advantage over Hartree-Fock methods 
that certain correlations are already build in. From 
the PESs we determine the activation barrier and pro- 
vide a new estimate of the rate constant utilizing pre- 
viously estimated quantities like the eg - electronic cou- 
pling Hda ~ D<eg\H\eg>A between donor D and ac- 
ceptor A molecular species, the outer-sphere contribu- 
tion to the reaction barrier eI^^, and the preequilib- 
rium constant K, 



equ • 



Although our excitation energy 
of Co(lll), AE^+\equ = 13,120cm-\ is in rather good 
agreement with the experimental value, we deemphasize 
the absolute numerical values of excitation energies. In- 
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lations exploiting that DFT usually stabilizes high-spin 
over low-spin states. If Co(II) has a low-spin ground 
state, this would also explain why the optical absorption 
^Eg ^Tig has not been observed. 

The paper is organized as follows. In section || we 
review further experimental and previous theoretical ef- 
forts in more detail. In section HI we introduce the DFT 



code SI ESTA and explain our computational methods 
(section HI A), as well as present the calculation of the 
PESs (section IIIB). In section [V we provide a new 



estimate of the spin- allowed hexaammine self-exchange 
reaction based on our insights from the PESs and dis- 
cuss our results. Finally, we summarize in section 



II. REVIEW OF PREVIOUS EFFORTS 

In this section further experimental facts about the 
single complexes and a more detailed review of past the- 
oretical efforts are presented. 

^3+ 




It is well accepted that Go{NH^) + [Co(III)] is a sta- ^IG. 1: Introduction of the ground state ^T^ [Mig] and first 
ble low-spin (S=0) compound, while it is assumed that excited state '^Eg pTig] of Co(II) [Co(III)], their irreducable 

representations, total spin S [fi], degeneracy g, and possible 
couplings (1) to (4). The energy ordering corresponds to the 
equilibrium Co-N bond lengths and rf of Co(II,S=3/2) and 
Co(III,S=0), respectively. 



Co{NH^)l+ [Co(II)] is high-spin (S=3/2). There are sev- 
eral reasons for this. First it is known from Co NMR 
studies that the related system, Co(i?2 0)|^, shows no 
sign of exchange with paramagnetic species fl^ and hence 
is low-spin. This is also expected to be true for Co(III), 
since it has an even larger ligand- field favoring low-spin. 
X-ray diffraction data of related crystals shows further a 
drastic difference of the Co-N bond distances (« 0.22A 
H, 0) between Co(II) and Co(III). From this it was 
concluded that the ligand-field of Co (II) must be much 
smaller resulting in high-spin. Additionally, the optical 
excitation spectrum of Co (III) could be fully character- 
ized including d — d spin- forbidden transitions. This is 
different than Co(II), where it could not be measured suc- 
cessfully. It is further experimentally supported that the 
Co^NH^)^^^^ self-exchange reaction occurs as outer- 
sphere. This is due to a rather slow ligand exchange rate 
of Co (III). There is considerable thermodynamic and ki- 
netic stability, which arises from effective cr-donation into 
empty eg-shellsp^. 

The theoretical effort is conveniently discussed in the 
context of the separable semi-classical transition-state 
model 



(2) 



In this equation, K^qu is the preequilibrium factor de- 
scribing the probability to form a precursor compound, 
Vef f is an effective nuclear attempt frequency to reach the 
transition state (TS), n^q < 1 is the electronic transmis- 
sion coefficient evaluated at the TS and averaged over all 
possible precursor compounds, r„ > 1 is the nuclear tun- 
neling factor, and the Boltzmann factor (classical Franck- 
Condon factor) gives the probability to reach the TS with 
activation free-energy {(3 = l/{kBT)). For the further 



states in figure ^ i.e. their group theoretical representa- 
tions, total spin S, degeneracy g, and possible couplings 

(1) to (4). The couplings (1) and (3) are due to spin-orbit 
coupling or thermal population at finite temperature T, 

(2) and (4) are mainly due to mixing of the electronic Cg- 
orbitals of Co(II) and Co(III) mediated by their ligands 
when a precursor state is formed. 

It has been known for about 50 years that extraordi- 
narily slow rates can often be attributed to small Franck- 
Condon nuclear overlaps. These originate from large 
changes in the metal-ligand bond length in the first coor- 
dination sphere Ip^. This seemed to apply also to the case 
at hand, if one assumes an adiabatic reaction [nei ~ 1) 
"neglecting" the spin-forbiddeness. However, a first de- 
tailed theoretical study by E. Buhks and co-workers us- 
ing perturbation theory in the weak spin-orbit coupling 
to admix spin-excited states shows that the electronic 
factor is very diabatic, k^i ~ 10^^ This makes the 
rate constant too small. 

Orgel[^ and by Stynes and Ibers[^ put out an 
idea that the reaction could involve thermally excited 
Co(II,S=l/2) making it spin-allowed. Larsson et al |2^ 
went further and argued that the excited states can be- 
come much lower in energy near TS. The important quan- 
tities are the energy difference between first excited state 
and ground state of Co(II) and Co(III) as a function of 
the Co-N bond length r 



AE^+( 



A m2-l-/„\ 



171/2 m \ 



(3) 
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In the case of Co(III) the lowering enhances some- 
what the admixture to the ground state, while in case 
of Co(II), "^Eg becomes even lower than ^Aig circum- 
venting the spin-barrier. Their argument is based on 
Born-Oppenheimer potential energy surfaces (PES) cal- 
culated both with ab initio Hartree-Fock (HF) and semi- 
empirical INDO-CI methods including e^ JT effects for 
both excited states. Unfortunately their result was not 
convincing, since the HF calculation predicts the wrong 
ground state for Co(III), while the INDO-CI excitation 
energy Ai?^+|eg„ = 2, 000 — 3, 000 cm^^ was much lower 
than the generally believed 9,000 cm^^[^, although 
there is no experimental evidence for this value. There 
has been no absorption found in this region. On the 
other hand, the excitation energy of Co (III) is known 
from experiment to be 13,700 cm~^[^ at the Co(III) 
equilibrium configuration. Since their calculated value 
AE^^lequ — 14,800 cm^^ is too large by an amount 
c « 1, 100 cm"\ the PES from INDO-CI were corrected 
according to 

AE^+ AE^+-c (5) 
AE^+ AE^++c. (6) 

It was noted that the semi-empirical INDO-CI method 
is generally very successful in calculating spectra of 
transition-metal complexes at fixed geometries, but not 
so in predicting molecular geometries [2^ . 

In 1991, Newton carried out new ab initio calcula- 
tions at both the SCF(UIIF) and correlated (UMP2) 
level. Using an empirical correction factor as large as 
c = — 6, 000 cm^^ (HF generally favors high-spin), he es- 
timated the excitation energies at the transition state to 
be AE^+ = 8,800 cni^i and AE^+ = 5,300 cm^^, and 
concluded that thermally excited pathways are not com- 
petitive to the spin-forbidden ground state pathway at 
room temperature (RT). This left the problem unsolved. 

III. AB INITO CALCULATIONS 
A. Method 

Our results of the inner-sphere activation barrier are 
based on PESs obtained from the fully ab initio code 
SIESTA||l|| based on density functional theory (DFT). 
SIESTA uses TrouUier-Martins norm-conserving pseudo 
potentials (25j in the Kleinman-Bylander form|g6j. For 
cobalt, we included spin-polarization and non-linear core 
corrections ||2^ to account for a spin-dependent exchange 
splitting and correlation effects between core and valence 
electrons, respectively. Relativistic effects are included 
for the core electrons in the usual scalar-relativistic ap- 
proximation (mass-velocity and Darwin terms) and by 
averaging over spin-orbit coupling terms, while no spin- 
orbit coupling is included for the 4s and 3d valence elec- 
trons. This has the computational advantage that spin 



of the total energy are less crucial when one is inter- 
ested in total energy differences. SIESTA uses a local 
basis set of pseudo atomic orbitals (PAO) of multiple 
type. The first-^ orbitals are produced by the method 
by Sankey and Niklewskip8[, while the higher-C orbitals 
are obtained from the split valence method well known 
from quantum chemistry. Polarization orbitals can also 
be included. We used a double-C basis set with polar- 
ization orbitals (DZP), as well as the generalized gra- 
dient approximation (GGA) in the version by Perdew, 
Burke and Ernzerhof |23| for the exchange-correlation en- 
ergy functional. The charge densities are calculated on a 
real space grid, where the fineness of the grid corresponds 
to an energy cut-off 80 Ry. 

Before going into the details of how we obtained the 
PESs plots, we outline the main ingredients of our calcu- 
lations. 

• Our calculations rest on the fixed-spin method 
within the spin-polarized DFT frame work. Two 
different Fermi energies, one for spin-up and one 
for spin-down, are adjusted in a self-consistent way 
in order to obtain the ground state of a desired to- 
tal spin. If a ground state of a certain spin is higher 
than the ground state of a different spin, then we 
know the former as the excited state and the lat- 
ter one as the true ground state. Nevertheless, 
there might be several states within a spin-manifold 
which differ in orbital symmetry. Since the DFT 
method is based on the variational principle, we 
obtain the state of a system with a certain total 
spin which has the lowest energy Our excita- 
tion energies are extracted from total energy differ- 
ences using the well-established self-consistent-field 
method (ASCF)|^. This gives generally reliable 
results for molecules since final state effects are in- 
cluded. 

• The ammines are treated as rigid bodies, i.e. the 
N-H bond lengths and the H-N-H bond angles 
stay the same throughout all our calculations. 
This is justified, since the vibrational frequencies 
of covalent bonds (N-H) are about one order of 
magnitude higher than the metal-ligand stretching 
frequencies 1^ and are considered average values. 

• We perform geometrical conjugate gradient opti- 
mizations (CGOs) of both Co(II) and Co(III) with 
the spin being fixed to a desired value. This is 
used to obtain the equilibrium geometries (EQ- 
GEOs) and the PES. The EQGEOs are optimized 
or relaxed geometries, where the forces on cobalt 
and all the ammines are below the chosen toler- 
ance O.OleV/A. As for the PESs calculations, we 
constrain the maximum size of a CG step to be 
O.OOlA in order to provide a large sequence of en- 
ergies. The resulting total energy can subsequently 
be plotted as a function of the Co-N bond length 
averaged over the six ligands. In the case of JT dis- 

j TVT 1 1 1 :_l 
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direction, rax, a-nd the average Co-N bond length 
in the equatorial plane, r^q, are more reasonable 
choices. 



B. Calculation of potential energy surfaces 

In the following we describe the procedure of how we 
obtained the PES plots in figures |^ and ^. 

Figure p| shows the ground state (S=0) and excited 
state (8=1) PESs of Co(III). For the ground state PES 
we started a CGO at the Co(II,S=3/2) EQGEO, de- 
noted as starting point (1), but with the spin fixed to 
S=0 and charge -|-3e. For the excited state PES (S=l) 
we started CGOs at both the EQGEOs of Co(III,S=0), 
starting point (2), and Co(II,S=3/2) (not shown for clar- 
ity). The vertical arrow indicates the optical excitation 
energy from the Co(III,S=0) EQGEO, where it is as- 
sumed that the excited state is initially not JT distorted, 
since the excitation is almost instantaneous. The dotted 
line in figure |^ corresponds to total energy versus the 
average Co-N bond length. Since Go(III,S=l) is JT un- 
stable, we also plot the total energy versus the average 
axial Co-N bond length, Vax, and the average equato- 
rial Co-N bond length, req, shown by dashed lines. One 
can see that essentially only r^q changes when starting a 
CGO from the Co(III,S=0) EQGEO. 

In order to obtain smooth PESs it is essential to treat 
the ammines as rigid bodies. If one relaxes all the 
atoms, the resulting PES would be much more compli- 
cated, i.e. the PES would display a sequence of short 
parabolic curves stemming from periods of contracting 
or stretching the N-H bond lengths, or from changing 
the H-N-H angles. The resulting PES would not resem- 
ble smooth parabola-like curves (at least for the non- 
JT distorted ground states) anymore, since the reaction 
coordinate would not simply be the change of the Co- 
N bond length. Although equilibrium energies of the 
constraint complexes might be higher than the all-atom- 
relaxed ones, this is not a problem, since we are interested 
in energy differences and finally in the sum of the Co(II) 
and Co(III) system energies (overall energy shift which 
does not affect the rate). 

Figure || shows the ground state S=3/2 and excited 
state S=l/2 PESs of Co(II) as a function of the Co- 
N bond length averaged over all six ligands. Using 
an analogue procedure, the Co(II,S=3/2) PES results 
from a CGO starting at the Co(III,S=0) EQGEO, (1), 
where we have fixed the spin at S=3/2 and the charge 
at +2e. The dotted line shows the Co(II,S=l/2) PES, 
which was obtained from starting at the non- JT distorted 
Co(III,S=0) and Co(II,S=3/2) EQGEOs, (2) and (3) re- 
spectively. They stay non-JT distorted. The vertical 
arrow indicates the optical excitation to the non-JT dis- 
torted Co(II,S=l/2) state. Starting geometry (4) is JT 
distorted and hence lower in energy than (1). Further re- 
laxation with the CGO method (dashed curve) lowers the 




Co-N distance [A] 



FIG. 2; PESs of Co(III) for spin S=0 (solid line) and S=l [h\ 
(dotted line) as a function of the average Co-N bond length. 
The numbers (1) and (2) indicate the starting points of the 
CGO. Since Co(III,S=l) becomes JT distorted, we also plot 
the PES along the axial, rax (dashed line), and equitorial, 
req (dashed-dotted line), average bond length. The vertical 
arrow indicates the optical excitation discussed in the text. 



(1) Co(II) 




Co-N distance [A] 

FIG. 3: PESs of Co(II) for spin S=3/2 [h] (solid line) and 
S=l/2 (dotted and dashed hues) as a function of the aver- 
age Co-N bond length. The numbers (1) to (4) indicate the 
starting points of the CGO. The dotted line shows the non- 
JT distorded PES, while the dashed line shows a JT distorted 
PES of Co(II,S=l/2). The vertical arrow indicates the optical 
excitation discussed in the text. 



non-JT relaxed energy (dotted curve). The non-JT dis- 
torted PES (dotted curve) is hence only meta stable. The 
energy of the Co(II,S=l/2) EQGEO is only 0.3 kcal/mol 
(less than fc^T) higher than the Co(II,S=3/2) EQGEO. 

Table | contains calculated equilibrium Co-N bond 
lengths of the ground and excited states and their ex- 
perimental analogous. For comparison it also shows the 
ones of the Co{H20)^^^^ ground states. As typical for 
the DFT method they are overestimated, about 3% for 
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TABLE I: Equilibrium Co-ligand bond lengths for cobalt hex- 
aammine and hexaaqua complexes in units of A. Jahn- Teller 
(JT) distorted values: rax denotes the Co-ligand bond length 
in axial direction, while req in the equatorial plane. 





calcd. [A] 


exptl. [A] 


Co{NH,)f^ 


JT distorted ^Tig 


2.03 
rax = 2.01 
r,g = 2.20 


1.96@ 


Co(7Vii-3)r 


"Tig 

'Eg 

JT distorted ^Eg 


2.22 
2.09-2.17 
rax = 2.42 
= 2.03 


2.16[|| 






2.02 


1.87H" 




'Tig 


2.25 


2.08-210ly]' 



''p.447 
''p.441 



Table compares theoretical and experimental optical 
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excitation energies. The excitation energy Eg 
for Co{NHz)^^ is not known from experiment but is 
very important. If this energy separation is low enough 
Co(II,S=l/2) is well-populated at ambient temperatures 
and the electron transfer reaction is spin-allowed. Al- 
though absolute excited state energies from DFT cannot 
be trusted, trends established from DFT are often cor- 
rect. From table || we can obtain the following. Since 
DFT underestimates optical gaps such as the Dq ligand- 
field, DFT favors high-spin compounds. For instance, the 
ligand-field splitting of Co(III) is 274 kJ/mol|35 while 
SIESTA gives 246 kJ/mol. In other words, excitations 
from low-spin to high-spin i^Tig <— ^Ag) are under- 
estimated while excitations from high-spin to low-spin 
{^Eg ^ '^Tig and ^Tig ^ ^Aig) are overestimated. Sim- 
ilar results were found for singlet-triplet gaps of phenyl- 
nitrene and other hypovalent systems |Q. 

Unfortunately, in the case of low-spin to high-spin exci- 
tations there is only one reliable value for Co(III) and no 
experimental value available for Co{H20)q^ in order to 
confirm the trend. Nevertheless, we can easily see that 
the calculated gap for Co{H20)f^ (7, 150cto~-^) has to 
be a lower bound analogous to Co(III). The argument is 
as follows. From ligand-field theory the excitation energy 

which is 
]. Taking 



Tig <— Ag is additive in the ligand field 
about 25% larger for NH3 compared to H20p^ 
a value Dq = 27AkJ/mol for Co(III) |3| and 200kJ/mol 
for Co{H20)q'^ we can correct the excitation energy 
of Co(III), 13, 700cm~^, by the Z3g-difference and obtain 
an approximate value for Co(iJ20)g^, 7, 500cto~^. This 
is clearly higher than the calculated value 7, 150cm^^. 
Having this trend established we can conclude that the 
excitation ^Eg ^ ^Tig of Co{NH3)l+, 3, 680cm-i, is an 
upper bound on the real value. 

A further confirmation of the quality of our excitation 
energy of Co(II) comes from the fact that our value is sim- 
ilar to Larsson's with INDO-CI (3, 100-4, lOOcm"! after 



TABLE II: Spin-forbidden transition energies in units of 



^Tlg ^ 'Ag 


Co{H20)l+ 


Co{NH3)l+ 


low-spin 


calcd. 


7.15 


13.23 


to 


exptl. 


(> 1.89[pl)" 


13.70^]'' 


high-spin 


Eg ^ llg 






high-spin 


calcd. 


8.46 


3.68 


to 


exptl. 


« 6.4|33l'= 




low-spin 


■^Tlg <— "Alg 


Fe{H20)r 






calcd. 


19.01 






exptl. 


I4.22Q'' 







"from Cobalt-59 NMR fraction of high-spin < 10"'*. This is not 
an order of magnitude estimate, but a true upper bound set by 
experimental resolution. 

''corresponds to the maximum of the optical absorption peak at 
T = 8K which stems from a vertical transition in line with the 
classical Franck-Condon principle. 

■^from a fit of a four parameter octahedral ligand-field theory (Dq, 
B, C, A) including spin-orbit coupling to a circular dichroism spec- 
trum of [Zn{H20)6]Se04] doped by Co{H20)l+ at T = 80K. 
Ligand-field theory is rather good for electronegative ligands such 
as oxygen. 

''from ions doped into beryl crystal 



TABLE III: Energy difference A£;^+'^+'(r) in units of 
10"* cm~^ at the equilibrium position r'^_^_f^2+) of Co(III,S=0) 
(Co(II,S=3/2)) and at the transition state r^ of the conven- 
tional spin-forbidden process. The values of r^ of this work, 
Larsson's, and Newton's are 2.09, 2.03, and 2.14 A, respec- 
tively. 







r"2+ 


rt 






rt 


this work" 


AE-'+ 


13,23 


9.09 


AE''+ 


3,68 


-1.19 


Larsson''|22] 




13,7 


12.0 




3,1-4,1 


^ 


Newton 




13,7 


8,8 




9,1 


5,3 



"DFT-GGA method, the AE^+ value is w.r.t. non-JT distorted 
Co{II, S=l/2). 

''INDO-CI method, empirically corrected energies by Eqs. 



with c = 1, 100cm ^, energy difli'erence AE^^{r^) is estimated from 
their Fig. 5 

_f IJHF-I-UMP2 method, empirically corrected energies by Eqs. 
(|,|) with c = -6, 000cm -1 



the JT distorted Co(II,S=l/2) equilibrium energy is only 
slightly higher than the non-JT distorted Co(II,S=3/2) 
equilibrium energy. Since the S=l/2 PESs are likely 
too high by a constant energy shift, Co(II,S=l/2) could 
be close to degenerate with Co(II,S=3/2) or even be 
the true ground state. Note that these considerations 
do not include multiplet-spitting due to spin-orbit cou- 
pling, as well as entropy effects on the energy (~ —TS = 
— kBTln{g), where g is the degeneracy of the multiplet). 



T71 InrJl \ 



In table III we compare the important quantity 
A£^^+/^+(r) from Eqs. for various bond lengths 

r between this, Larson'sp^ and Newton's||^ work. 
^^3+/2+^y.o^^^^-j g^j.g excitation energies out of 

the equilibrium states Co(III,S=0)/Co(II,S=3/2) to 



/TTT n 
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energy difference at the spin-forbidden transition state. 
The excited state Co(III, S=l) is always much too 
high in energy to be relevant and is not important 
for further consideration. However, the value and sign 
(!) of Ai?^+(rt) decides whether the process will be 
spin-forbidden (allowed by weak spin-orbit coupling) 
(A£;2+ > 0), or whether it is spin-allowed (Ai?^+ < 0). 
The former process can be described by second order per- 
turbation theory |6[ |, using Co(II,S=l/2) as an vir- 
tually excited state coupled to at the spin- forbidden TS, 
the latter one is through direct coupling and does not 
involve spin-orbit coupling ||2^. As one can see our value 
is Aii'^^ (r'l' ) < and hence the reaction along the lowest 
energy pathway is spin-allowed. 

Figure ^ shows the PES of the total 2 Co-system, i.e. 
the sum of the single complex energies. The minimum 
in the left top corner at (r^ — 2.03A,rj^ = 2.22A) corre- 
sponds to equilibrium CoL{III)Cofj{ir), where the sub- 
scripts L and R stand for "left" and "right" , respectively, 
and simply distinguish the two complexes. The minimum 
in the right bottom corner at {vl = 2.22A, rR = 2.Q3A) 
corresponds to equihbrium Col{I I)C or{I 1 1) . The vari- 
ables tl and are the Co-N bond lengths of the left and 
right complex. For simplicity we can only show the non- 
JT distorted complexes, since a single JT distorded com- 
plex depends on two variables, rax and r^q. This would 
resolve in a total energy which depends on more than 
two variables and cannot be plotted. Part A is plotted 
using the ground state Co(III,S=0) (solid line in Fig. |^ 
and Co(II,S=l/2) for Co(II)-N bond length ra < 2.125A 
(dotted line in Fig. |) and Co(II,S==3/2) for > 2.125A 
(solid line in Fig. ||). This is an example (by not includ- 
ing JT distorted complexes of lower energy) of a spin- 
allowed reaction. However, if we restricted ourselves to 
the Oh group, this would be the lowest energy pathway. 
Such a reaction starts out at equilibrium Col{IItS = 
3/2)Coii{III, S — 0), then activated by thermal fluctu- 
ations changes to the Col{II, S = 1/2)Cor{III, S = 0) 
PES at {rL = 2.1251, rj? = 2.031) through spin-orbit 
coupling in only first order (presumably adiabatic) . Com- 
ing back to figure |, part B shows the PES of Co(III,S=0) 
and Co(II,S=3/2) and hence describes the conventional 
spin-forbidden reaction with an energetically higher tran- 
sition state. 

In the following wc have to obtain the activation bar- 
rier of the reaction. In figure || the transition states 
for spin-allowed (lower graph) and spin-forbidden (up- 
per graph) are plotted. The graphs are the cross-section 
of the 2Co-PES from figures ^ A and B along the diago- 
nal r ■=rL= TR, i.e. E{rL = r,rR = r) = Eco{iii){r) + 
Eco{ii){r), where the Co-N bond lengths and tr are 
the bond lengths of the "left" and "right" complex. The 
transition state is defined as the saddle point. Using 
non-JT distorted Co(II,S=l/2) leads to a lowering of the 
spin-forbidden transition state by O.SeV = 6.92kcal/mol. 
The spin-allowed and spin-forbidden TSs are at r = 2.05 
and 2.ll, respectively. In the next section we use the 
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FIG. 4: Contour plot of the PES of total 2 Co-system with 
contour level spacing 0.03 eV = 0.69 kcal/mol. Only results 
from the non-JT distorted complexes are shown. A) Example 
of a spin-allowed reaction using Co(III, S=0) and non-JT dis- 
torted Co(II, S=l/2) for Co(II)-N bond lengths r2 < 2.125A 
and Co(II, S=3/2) for ra > 2.125A (see Fig. |). Within 
the Oh group this is the lowest energy reaction pathway. B) 
PES of conventional spin-forbidden reaction originating from 
Co(III, S=0) and Co(II, S=3/2). 



mine self-exchange reaction. 

Finally a few things are important to keep in mind. 
Application of the DFT method to transition metal ions 
is tricky, even more when one is interested in excited 
states. There are several well-known deficiencies one has 
to consider. First, it is well known that DFT suffers from 
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E [eV/precursor complex] 

-5600.2 - 



E (V, S=l/2) 



(2) 



e"(V, S=1/2) 




2.15 2.2 

rL=i-R [A] 



(3) 



E°(V, S=3/2)= -5601.39 



FIG. 5: Cross-section of 2Co-PES, written as E(total charge 
in roman letters, total spin), along tl = tr,. The TS is defined 
as the saddle point. Using Co(III,S=0) and non-JT distorded 
Co(II,S=l/2) lowers the TS by O.SeV/precursor complex = 
Q.Q2kcal/mol compared to using Co(III) and Co(II,S=3/2) 
(arrow (1)). Arrows (2) and (3) indicate the two activation 
barriers being the difference between the energy at the tran- 
sition state (superscript f) and at the EQGEO (superscript 
0). 



concerns mainly the localized d-orbitals of cobalt, where 
the large d—d Coulomb interaction introduces local corre- 
lations that are not captured properly by the GGA func- 
tional. The second deficiency concerns optical gaps which 
are generally underestimated |39[ |40). In our case the 
excited states involve the ligand-field spitting. This, how- 
ever, is turned to our advantage by using it to deduce a 
trend. It is utilized in the next section to obtain an essen- 
tial lower bound on the rate constant. More importantly, 
since DFT is a ground state theory, we expect that ex- 
citation energies obtained from total energy differences 
(ASCF method) are rather reliable. Very good results 
have been obtained for optical spin (singlet-triplet) Q 
and charge]^ excitations. Furthermore, total energies 
are expected to be better for strongly cr-donating NH^ 
than for weakly 7r-donating H^O ligands, because the 
NH^, ligands are less electro-negative and bind more co- 
valently with cobalt Cg orbitals. Covalent molecular-type 
systems are well described by the DFT method. Besides, 
covalency screens the on-site repulsion of the localized 



IV. RESULTS 

In this section we give an estimate for the hexaam- 
mine self-exchange rate using previously estimated elec- 
tron transfer parameters and results from last section. 
For now, we take the energies obtained from the DFT 
calculations literally. For instance, we assume that high- 
spin Co(II,S=3/2) is the groundstate according to our 
calculation (see Fig. ||). However, since DFT stabilizes 
high-spin Co(II,S=3/2) over low-spin Co(II,S=l/2) ac- 
cording to our trend, this rate estimate is a lower bound. 
This is because we first have to thermally excite to the re- 
action intermediate, Co(II,S=l/2), which lowers the rate. 
With regard to our trend, Co(II,S=l/2) could be the true 
groundstate. 

The initial rate constant to the intermediate state 
Co(II,S=l/2)Co(III,S=0) is denoted by h. In ther- 
mal equilibrium the return reaction with rate constant 
k-i back to the groundstate Co(II,S=3/2)Co(III,S=0) is 
equal to ki . From the intermediate state the spin-allowed 
electron transfer can occur with rate constant ket 



Col{II, S = i/2)CoR{III, S^O) 

Codll, S = 1/2)Cor{III, S = 0) 
CoL{ni, S ^ 0)Cor{II, S=l/2) 



ki 



(7) 



The resulting total rate constant is k. Since the two steps 
from Eq. |^ are incoherent due to relaxization along the 
surfaces, the total rate k can equivalently be described as 
starting out at the thermally populated (with probabil- 
ity P) intermediate, Col{II,S = 1/2)Cor{III, S = 0), 
from were the spin-allowed electron transfer can occur to 
Col{III,S = 0)Cor{II,S = 1/2) with rate ket- This 
does not require rate constant ki. Hence, the total spin- 
allowed reaction rate constant is 



k = Pket ^ 



9 e 



(8) 



with ket being the electron transfer rate Eq. (||). 
The prefactor or probability P applies in the case that 
Co(II,S=l/2) is the excited state, which has first to be 
thermally populated. AG* = AE* - TAS is the free 
energy of excitation and Z is the partition function. 

There are two main differences from previous rate 
estimates: 

• Co(II,S=l/2) and Co(II,S=3/2) are nearly degen- 
erate, i.e Co(II,S=l/2) is shghtly higher in energy 
than Co(II,S=3/2) but is probably overestimated. 
In order to get a lower bound of the rate we hy- 
pothetically trust the DFT energies and calculate 
AG* = 0.3kcal/mol + kBTln{12/4) = 3.1kcal/mol 
from the energy and the entropy differences be- 
tween Co(II,S=l/2) and Co(II,S=3/2). Hence, the 
prefactor of equation (ph is 0.09 using the degen- 
eracies g from figure nPand neglecting multiplet- 
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• The TS of the Co(II,S=l/2)Co(III,S=0) system 
is lower than the TS of Co(II,S=3/2)Co(III,S=0). 
The inner-sphere contribution to the activation 
barrier, the energy difference of the 2 Co-system be- 
tween the TS and equihbrium using Co(II,S=l/2), 
is only Ej^ = 0A2eV = 9mkcal/mol. For the 
equihbrium energy of Co(II,S=l/2), we take the JT 
distorted value, for TS we use the non-JT distorted 
one. This seems reasonable, because the acceptor, 
Co(III,S=0), is non-JT distorted. In order to have 
sufficient nuclear overlap, the donor should be non- 
JT distorted at TS, too. 

The other parameters are only slightly modified. This 
is mainly due to using a different nuclear frequency for 
Co(II), i.e. MEg) for Co(II,S=l/2) instead of ^Aig) 
for Co(II,S=3/2): 

The product Keq„i^e// = 1.4 - 3 • 10"(Ms)-i has 
been estimated before for the spin-forbidden process by 
Sutin [|j and Newton|^, |l^ using a preequilibrium con- 
stant Kequ — 0.013A/^^ind an effective nuclear fre- 
quency Veff = 347cm~^||^. Since the case at hand is 
sHghtly different (low-spin Co(II, 8=1/2)), we redeter- 
mine i^eff- The average harmonic frequency is given 
by|, H, ID 



2 

solv solv 



.El 



E, 



solv 



El 



(9) 



where 



y| is the reduced inner- 



respectively we obtain Vin = 320cr7i ^ = 
For the solvent, we use a typical value 

irl 



sphere nuclear frequency of Co(II) and Co(III). Using the 
experimental values V2(,Eg) = 255cm~^ and 1/3(^13) = 
494cm- i|l7 
9.6 • IQi^s- 

30CTO-1 = 0.9 • 10i2s-^[^. We do not attempt to es- 
timate i>2/3 from the PESs due to complications from 
JT distortions in the case of Co(II,S=l/2). The outer- 
sphere, = 28 /4kcal/mol = Ikcal/mol, was previ- 
ously determined by Buhks and colleges applying the 
Marcus-Levich continuum model for the solvent |l^, p^ . 
We neglect entropy contributions, which are quite small 
for water near room temperature and vanish for the 
inner-sphere system, if the (harmonic) vibrations are the 
same for the activated complex and the reactants[^. Us- 
ing these parameters, one obtains v^ff — 245cm-^ — 
7.3 • IQi^s-i and K^quVeSf = 9.5 • W°{Ms)-^. 

The electronic transmission coefficient for the Cg- 
transfer was determined by Larsson and colleges to be 
weakly adiabatic, k^i ~ 0.5[p2[, using parameters for 
the spin-forbidden process and standard Landau-Zener 
theory Extended Hiickel theory was applied to 

calculate the donor(D)-acceptor(A) electronic coupling, 
Hda = D<eg\H\eg> A, for different precursor complexes 
(apex-to-apex, apex-to-edge, apex-to-side). A subse- 
quent orientational averaging gave Hda ~ 200cm-^. 
Their statistical analyses assumes that each complex can 
rotate independently and that each configuration covers 



^"i/e// and Eqs. (7- 
we redetermine the electronic trans- 



lA (van der Waals contact between first solvation shells), 
the Co-N distances of both complexes were fixed at 2. 06 A 
reasonably close to our TS value, 2.05A. A similar pro- 
cedure was applied by Newton Utilizing Hda from 
above, our parameters e\^i^, E^ 
11) fromRef.fl, _ 
mission and nuclear tunneling factors and obtain adia- 
baticity k^i = 0.73 and weak nuclear tunneling r„ = 1.9, 
respectively. 

Finally, the activation barrier has inner-sphere and 
outer-sphere contributions and is lowered by the average 
electronic coupHng at TS, i.e. w = eI^ + eI^^^ - 
Hda, and is given by = 0.70ey = 16.Qkcal/mol. 
This leads to a Boltzmann factor e^^'-^ = 2.1 • 10"^^. 
Including all the calculated charge transfer parameters 
in Eq. (0) gives a rate constant of 6 • 10~'^(Ms)-^, which 
is about 2 orders of magnitude larger than experiment. 
Possible sources of errors are discussed in section ^ 

In the following, we again want to stress the main dif- 
ferences between our treatment of a spin-allowed process 
and the conventional spin-forbidden reaction by Buhks 
et al. and Newton|^. The main difference is that in 
our case the low-spin Co(II,S=l/2) complex has a lower 
energy than the high-spin Co(II,S=3/2) near the TS as 
opposed to Buhks et al. and Newton. In their case, 
the spin excited state Co(II,S=l/2) is much higher in 
energy and is only virtually coupled to by weak spin- 
orbit coupling. This leads to extremely small rates of 
order 10-^''(Afs)^^. In table [V we show for further il- 



lustration the individual parameters used to calculate the 
rate according to Eqs. (|^) and ^ and compare to New- 
ton's spin-forbidden and hence diabatic groundstate reac- 
tion and thermally excited adiabatic pathway. The spin- 
forbidden reaction is significantly lowered by the weak 
spin-orbit coupling (ke/ « 1), while the excited alterna- 
tive is thermally not accessible and leads to even smaller 



rates. The ^E„ 



Tig excitation energy is as large as 



9, 100cm , while in our case it is only 3, 680cm 



V. DISCUSSION 

At this point, there are two main possibilities why 
there is a 2 order of magnitude disagreement between 
our, 6 • 10~"^(Afs)-^, and the experimental rate, 5 • 
lO-^(Afs)-^. Either something is not correct with the 
theoretical estimate, or the experimental rate constant 
is too small. In the more likely case that the theory 
misses some details, the sources of possible errors are the 
^Eg <~ "^Tig excitation energy and hence the prefactor 
of Eq. (||). Increasing the energy of the Co(II,S=l/2) 
PESs w.r.t. Co(II,S=3/2) by a constant energy shift 
2, 560cm- ^ and hence increasing the excitation energy 
to 6,240cm-^, produces the experimental value. This, 
however, is in conflict with the trend of excitation ener- 
gies established from the DFT method. 
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TABLE IV: Rate parameters (see Eqs. (g) and (g)) used in this (spin- allowed) and Newton's(spin-forbidden and thermally 
excited) work. The rates are compared with experimental estimates. (!) indicates main differences between Newton's and our 
work. 





g/z 


K 






r„ 




rate fc[(Afs)"^] 


exp.[3, 1 1^ 














10" ' - 10"" 


this work 


0.28 


0.013 


245 


0.73 


1.9 


16.3" 


6 ■ 10"^ 


Newton"!]: 
















thermally excited 




0.013 


245 


0.67 


2.4 


30.7(1)'= 


2 ■ 10"" 


spin-forbidden 




0.013 


347 


10~*(!)'' 


9 


24.4 


4-10-^° 



"includes AG* from Eq. (|8|) 

"individual parameters are presented to the best of our knowledge 
■^due to a large ^Eg <— ^Tig excitation energy of 9, lOOcm"^ 
''due to a spin-orbit reduction factor of 7^ = 1.8 • 10~*. 



originate from a more complicated cross-over from the 
JT distorted Co(II,S=l/2) equilibrium complex to ther- 
mally excited non-JT distorted Co(II,S=l/2) near TS. 
The cross-over could be rather unlikely, since the sym- 
metry changes. 

On the other hand, possible experimental issues are 
beyond our realm of knowledge. Although a difhcult ex- 
periment, it has been thoroughly studied over decades. 
Nevertheless, generally forgotten, neglected or underes- 
timated side-reactions lead to an underestimation of the 
rate constant. In particular, the existence of high-spin 
Co(II,S=3/2) is challenged by our analysis. This leads to 
the question, if the assumption of high-spin Co(II,S=3/2) 
was wrongly made when deriving certain rates of side- 
reactions. Having said this, an error of about 2 orders of 
magnitude in the rate estimate is not as bad as it may 
sound. Due to the activated nature of the reaction, the 
exponential dependence of the rate on the activation bar- 
rier and excitation energy makes it very sensitive to small 
errors. 



VI. CONCLUSION 

Within Marcus theory of charge transfer the rate 
constant of the hexaamminecobalt-self-exchange reac- 
tion was redetermined. We utilized the DFT code 
SIESTA to calculate Born-Oppenheimer potential en- 
ergy surfaces and spin-excitation energies and use pre- 
viously determined parameters. The main differences 
from former work is the near degeneracy of Co(II,S=3/2) 
and Co(II,S=l/2). Furthermore, we observed a dras- 
tic lowering of the activation barrier (~ 6.9kcal/mol) 



for the reaction pathway involving and Co(III,S=0) and 
Co(II,S=l/2). This led to a rate constant of order 
6 • 10~'^(Ms)"^ which is 2 orders of magnitude faster 
than experiment. Possible sources of errors are out- 
lined in section ^ and having most likely to do with 
the neglect of the proper cross-over treatment from 
Jahn- Teller distorted equilibrium Co(II,S=l/2) to non- 
Jahn- Teller distorted Co(II,S=l/2) near the transition 
state. The good quality of our energetics involved in 
the rate constant evaluation is indicated, first because 
our calculated excitation energy 13, 230cto"^ of Co(III), 
^Tig <— ^Ag, agrees well with the experimental value 
13,700cm~^. Second, our excitation energy 3, 680cto~^ 
of Co(II), ^Eg <— ^Tig, is close to Larsson's value from 
INDO-CI (3, 100 - 4, lOOcm-i). This energy is in par- 
ticular important for charge transfer. Unfortunately, 
the corresponding optical absorption peak could not be 
found experimentally. This may simply be explained 
by the spin-forbiddeness and vibrational broadening be- 
cause of differences in equilibrium Co-N bond distances 
and a Jahn- Teller unstable excited state. On the other 
hand, our analysis questions the existence of the high- 
spin ground state Co(II,S=3/2) and hence can provide 
an alternative explanation for the absence the the ab- 
sorption peak. 
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